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Abstract. The aim of this paper is to investigate the conditions for the coexistence of phases 
in a Lennard Jones fluid. The calculation has been performed within the virial developement 
method, and as a result, a simple approximate relation has been obtained between the number 
densities of two coexisting phases and the interparticle interaction potentials in them. The 
results of this work could have implications for modelling of giant planet interiors. This has 
become important, due to the discovery of more than 300 extrasolar planets. 



1. INTRODUCTION 

In the last 5-6 years of the XX century,interest in planetology has drastically risen. 
This is due partially to results obtained within our planetary system, but even more 
so to discoveries of extrasolar planets. Until the end of November 2008., according 
to data at http://exoplanet.eu , 329 exoplanets have been detected and 269 of them 
have masses M < 5Mj where Mj is the mass of Jupiter. This testifies about the 
interest of modelling the internal structure of the giant planets. For a recent study 
see, for example, (Vorberger et al.,2007 or Nettlemann et al.,2008). 

The aim of this contribution is to present preliminary results on the conditions 
under which two phases in a Lennard- Jones fluid can be in equilibrium. Phases are 
defined as regions of the parameter space within which properties of a material are 
uniform. The condition for their equilibrium is the equality of pressures, temperatures 
and chemical potentials. In this contribution only the equality of pressures and tem- 
peratures will be considered; accordingly the results will be only preliminary. The 
equality of chemical potentials will be included in future work. 

2. METHOD 

The virial development of the equation of state is a method in which the equa- 
tion of state (EOS) of a fluid can be expressed as a power series in the density,and 
the coefficients take into account the interactions present in the system, in which an 
increasing number of particles takes part (e.g. Reichl 1988). 

^Presented at the 15*'' National Conf.of Astronomers of Serbia, Belgrade, October 2-5. ,2008. The 
proceedings will appear in Proc. Astron. Obs. Belgrade. Disregard the volume and page numbers. 
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The mathematical form of the EOS of a fluid within the virial development is 

All the symbols on the left side of Eq.(l) have their standard meanings, while on the 
right hand side, a; are the so called virial coefficients, A is the thermal wavelength and 
V is the inverse number density of the system v = V/N. The thermal wavelength is 
given by (for example Reichl,1988) 

where h is Planck's constant and m the particle mass. The LJ model potential has 
the form 

u{r) = 4e [{^r - i^r] (3) 

The symbol e denotes the depth of the potential,while a is the diameter of the molec- 
ular "hard core". The first coefficient in Eq.(l) is a; = 1 , while the second one is 
given by 



0-2 



(T) = -(2/3)7rArA/3 exp(-/3u)(— )r3dr (4) 

(Maitland,G,C.,Rigby,M., Smith. B.E. and Wakeham,W.A.,1987). In this and other ex- 
pressions P = I/UbT ,T is the temperature and fc^ is the Boltzmann constant.lt can 
be shown that the second virial coefficient for the Lennard- Jones potential is given by 

oo 

a2{T*) = boJ2^j{l/T*Y'^+'^/^ (5) 
j=o 

where 

T* - ^ (7) 

and 6o = (27r/3)A^yi(T'', where Na is Avogadro's number. The chemical potential of a 
fluid is given by (Hill,1987): 

n 

In(nA^) -I- / d-f / dr4nr'^u(r)g(r) (8) 

ksT Jo Jo 



where n is the particle number density and g{r) is the radial distribution function. The 
problem with the calculation of ji is the determination of 5 (r), which is a complicated 
task in statistical mechanics. In the future, it will be attempted to insert the form of 
g{r) proposed in (Morsali et.al.,2005) into Eq.(8). 
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3. THE CALCULATION 

In order to render the EOS within the virial development physically applicable, 
the power series in Eq.(l) has to be convergent. This means that the results are 
applicable only under the condition 

— =nX^^l (9) 

V 

Inserting Eq.(2) into Eq.(9),it follows that the virial development is applicable 
under the conditions 



n -< 



/ m \ 3/2 




\2^) 





(10) 



Denote the two mutually non-interacting phases which make up the system by " 1" 
and "2". Applying Eq.(l), gives 

°° \3 

p,=n,kBT,Y,MT){-y-' (11) 
1=1 

and 

°° \3 

P2 = n2kBT2y^Cr{T){^Y-' (12) 

Inserting the conditions of the equality of pressures and temperatures needed for the 
coexistence in equilibrium of the two phases ,one gets the following expression for the 
ratio of densities in them 

Taking into account that the equality of the temperatures implies the equality of the 
thermal wavelengths, this expression can be transformed into the following form 

oo oc 

J2ai{T)n[X'-' = ^c,(r)n^A'-i (14) 

1=1 r=l 

Limiting the sums to the first two terms, it follows that 

aiui + a2n^A = Cin2 + 02^2'^ (15) 

Taking into account that the first virial coefficient is l,and introducing m — 722 = x, 
the last expression can be solved to give 



X = (l/2Aa2) -1 - 2Aa2n2 + y/l + 4a2n2A(l + C2n2A) (16) 

This result is mathematically simple, but physically interesting. It gives the difference 
between the number densities of two phases coexisting in equilibrium,expressed in 
terms of the density of one of the phases ,the thermal wavelength and the second 
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virial coefficient in both of them. Expression (16) can further be transformed to give 
finally 

1 



1 - V'rT4^^^^^^A(TT^^^^^ (17) 
zAa2 L J 

which is positive for 4o2n2(l + 02/12 A) >- 0. 

This result is even more physically interesting. It represents a link ( of second 
order, because only the second virial coefficients are taken into account) of the number 
densities in two phases of a LJ fluid, expressed as a function of the temperature and 
the interaction potentials in both of them. This result is general,in the sense that 
virial coefficients for any potential can be inserted in it. 

Turning to the case of a i — J fluid,the values of the first few coefficients 7j are: 
7o = 1.733,71 = —2.564,72 ~ —3.466. This means that the explicite expression for 
the second virial coefficient is 

27r 



a2 = -f^ACT^ 



l.T33(— )'/'-2.564(— )»/'-3.466(— )VV 



(18) 



Inserting Eq.(18) into Eq.(17),one would obtain an expression "linking" the number 
densities of the two phases coexisting in equilibrium with the parameters of interpar- 
ticle potentials in them. How could this result be applied in studies of the interiors of 
the Jovian planets? Within any celestial body, the particle number density increases 
with increasing depth. Therefore, the parameters of the L — J potential would also 
be density dependent, which implies that the virial coefficients would also be depth 
(an density) dependent. This means that the possibility of phase coexistence would be 
density dependent. A closely related problem is the density dependence of the phase 
transition pressure in a system with a L — J potential. These and other issues related 
to problems of the behaviour of a L — J fluid under high density will be studied in 
future work. 
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